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Abstract. The angular analysis of exclusive rare B-meson decays via intermediate vector 
mesons V into 4-body final states of two pseudo-scalars Pi , P2 and a pair of light leptons 
L = e, fi offers a large set of observables. They can be used to test the electroweak short- 
distance couplings in the Standard Model and to search for New Physics. The two kinematic 
regions of low and high dilepton mass depend on short-distance physics in complementary 
ways and can be expanded in powers of Aqcd/ti,;,. These expansions guide towards suitable 
combinations of observables allowing to i) reduce the hadronic uncertainties in the extraction 
of the short-distance couplings or ii) test the lattice QCD B — > V form factors in short-distance 
independent combinations. Several such possibilities of CP-averaged and CP-asymmetric (T- 
even and T-odd) quantities are presented for B^ — > K*°{-^ K~n + ) + it and time-integrated 
CP-asymmetries without tagging for B S ,B S — > <f>{— > K~ K + ) + ii decays in view of the latest 
B-factory and CDF results and the forthcoming LHCb measurements. 



1. Introduction 

The b — > s+ll induced flavour changing neutral current (FCNC) decays of the B- meson are well 
known to provide sensitive probes of the electroweak short-distance couplings of the Standard 
Model (SM) and to test scenarios beyond (BSM). Hitherto existing experimental results from 
the .B-factory experiments Babar [1] and Belle [2] as well as the Tevatron experiment CDF [3] 
are in the ballpark of the SM predictions and provide constraints on the couplings Cg^o related 
to 4-Fermi operators (s7u-Pl6) (^7^{1, 75}^) [4, 5]. Further experimental progress is expected 
soon from the final analysis of the full data sets of Babar and Belle, as well as the CDF one 
with at least doubled statistics. Eventually, LHCb will dominate statistically for final states 
containing charged particles only, already with a rather low luminosity of a few fb" 1 [6]. The 
Super- £? factories are expected to contribute to I = e, just as the l?-factories do presently. 

Especially the angular analysis of the 4-body final state of B — > V(^- P\P-i)-\-lt decays offers 
a large number of observables in the fully differential distribution [7, 8] . Here the intermediate 
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Figure 1. Kinematic variables of 
£0 ^*o^ K--K+) + U decays: 

i) the (^)-invariant mass squared q 2 , 

ii) the angle 9i between £ = t~ and B 
in the (££) center of mass (cm.), Hi) 
the angle 9k* between K~ and B in 
the (K~tt + ) cm. and iv) the angle <j) 
between the two decay planes spanned 
by the 3-momenta of the (Kir)- and 
(££)-systems, respectively. 

V is assumed to be on-shell in the narrow-resonance approximation which restricts the number 
of kinematic variables to four 4 . Using B^ — > K*°(-^- K~ir + ) + U for illustration, they might be 
chosen as depicted in figure 1. 

The differential decay rate, after summing over lepton spins, factorises into 

~? Til a 12 = J i sin2 °K* + J\ cos 2 9 K * + ( J| sin2 °K* + J% cos 2 0* . ) cos 26> £ 

3 dq z a cos 6^ a cos 6^* d(p 

+ J3 sin 2 Ok* sin 2 cos 2cj) + J4 sin 26^* sin 20i cos + J5 sin 20k* sin (9^ cos 
+ ( Jg sin 2 0ft-* + Jg cos 2 Ok* ) cos + J7 sin 2(9^* sin 9i sin 

+ Js sin 29k* sin 2$^ sin + Jg sin 2 Ok* sin 2 sin 2(f), (1) 

that is, into g 2 -dependent observables 5 J\ (q 2 ) and the dependence on the angles 9?, Ok* and 
4>. No additional angular dependencies can be induced by any extension of the SM operator 
basis [11] as found by [12, 13]. The following simplifications arise in the limit mi — > 0: Jf = 3J|, 
Jf = - J| and = 0. 

The differential decay rate d A T of the CP-conjugated decay B® — > K°*(— >■ iT + 7r~) + is 
obtained through the following replacements 

J l,2,3,4,7 ^1,2,3,4,7^ ^ "M. J 5,6,8,9 ^ ~ ^5,6,8,9 ~> ~M> ( 2 ) 

due to £ £ ^ 9i ^ 9i — tt and — > — 0. The CP-violating (weak) phases 5^ are conjugated. 

The angular distribution provides twice as many observables (J| and J|) when the decay 
and its CP-conjugate decay are measured separately. This doubles again if the £ = e and fi 
lepton flavours are not averaged. Notably, CP-asymmetries can be measured in an untagged 
sample of l?-mesons due to the presence of CP-odd observables (i = 5,6,8,9) [7]. Moreover, 
T-odd observables ~ cos 5 S sin 5w (i = 7,8,9) are especially sensitive to weak BSM phases 5w 
[10, 14] contrary to T-even ones ~ sin<5 s sin 5w (i = 1, . . . ,6), since the CP-conserved (strong) 
phase 5 S is often predicted to be small. Note, that in the SM CP-violating effects in b — > s are 
doubly-suppressed by the Cabibbo angle as lm[V ub V* s / (y tb V t *)] « 77A ~ 10~ 2 . 

4 The off-resonance case has been studied in [9]. 

5 Possibilities to extract g 2 -integrated J\ from single-differential distributions in $£, 6 k* or (f> can be found in [10]. 




The observables j{ are bilinear in the transversity amplitudes Aa ,R (a = 0, _L, || in the limit 
mi — > 6 ) [15]. Here L, R refer to the chirality of the lepton current. The kinematic dependence 
on the lepton mass is suppressed by mj/q 2 . Hence, non- negligible effects from finite lepton 
masses at q 2 > 1 GeV 2 arise only in BSM scenarios, as for example in the MSSM at large tan /3 
due to neutral Higgs boson exchange, see for instance [16]. 

Depending on q 2 , present theoretical predictions suffer from long-distance dominated (qq)- 
resonance background (q = u, d, s, c) induced by current-current- and QCD-penguin-operators of 
the AB = 1 effective Hamiltonian. Compared to the leading c-current-current contribution, the 
■u-current-current one is suppressed by small CKM elements and the QCD-penguin contribution 
by small Wilson coefficients. The light (u, d, s)-resonances affect the FCNC b — > s + 11 decay at 
very low q 2 < 1 GeV 2 . The narrow (cc) -resonances J/ip and ip' are vetoed in the experimental 
analysis, while the broader (cc)-resonances appear at q 2 > 14 GeV 2 . 

This singles out a low-g 2 and a high-g 2 window which are usually chosen in the intervals 
q 2 G [1, 6] GeV 2 and q 2 > 14 GeV 2 , respectively. They correspond to large and low recoil of the 
K* -meson, as indicated by the scaling of its energy in the U-meson rest frame: Ex* ~ nib/2 
versus Ek* ~ rriK* + Aqcd- The heavy quark limit is combined at large recoil with the large 
energy limit as an application of QCD factorisation (QCDF) [17, 18] whereas at low recoil an 
operator product expansion (OPE) can be performed [19, 20] in combination with HQET form 
factor relations [21]. Revealing the symmetries of QCD dynamics in both regions, the expansions 
reduce the number of hadronic matrix elements. Many works have been inspired in the past 
decade in order to exploit these symmetries to reduce the hadronic uncertainties in the exclusive 
decays. Particularly, the angular analysis of the 4-body decay B — > U(— > P1P2) + ££ proved 
most fruitful. For comparison, B — > P + li decays offer only two additional observables beyond 
the branching ratio, the lepton forward-backward asymmetry ^4p B and the flat term Ffj [16] . 

2. B° -»• K°*(^ K-tt+)+££ 

The experimental situation in regard to lepton final states is twofold: On the one hand, the 
(Super-) .B-factories are able to measure both £ = e and {i modes separately. On the other 
hand, most detectors at hadronic collider machines prefer £ = fi, however, LHCb might be able 
to study £ = e, too [22]. 

2.1. Low-q 2 / Large recoil region 

Using QCDF in the low-g 2 region, the seven B — > V QCD form factors (V, A),i,2, ^1,2,3) reduce 
to two universal form factors £mi [23]. Furthermore, the B — > K* + it amplitudes factorise 
[17, 18] 7 , where the numerically leading contributions to the transversity amplitudes are [15] 

^occj^xfii, (3) 

with short-distance coefficients c j^'jp The transversity amplitudes have been used to construct 
the observables [15] and A^'^ [26, 27] with reduced hadronic uncertainties and an 

6 Two more amplitudes contribute for mi 7^ 0: At (time-like) and in the presence of scalar operators, As [12]. 

7 One might use equivalently soft collinear effective theory (SCET) [24, 25]. 



Table 1. The SM predictions of the g 2 -integrated CP-asymmetries of — > K°*(-^ K~ir + )+££ 
decays (q 2 £ [1,6] GeV 2 ) and their relative uncertainties from the form factors £_mi and the 
renormalisation scale ^ [10] (see also [12]). Also shown are the ranges of the CP asymmetries 
after applying the experimental constraints at 90 % C.L. for the generic model-independent 
BSM scenario [10]. Note that the SM predictions are given in units of 10~ 3 . 
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improved sensitivity to the chirality flipped operators O' 7 s l0 . Instead of fitting the observables 
Jf from the angular distribution, the authors of reference [26, 27] pursue the possibility to 
fit directly the transversity amplitudes which necessitates the identification of their rephasing 
properties. Moreover, the authors investigate the experimental reconstruction uncertainties of 
their method at LHCb, including among the theoretical uncertainties also lacking subleading 
QCDF contributions by simple power counting. 

The original analysis of the CP asymmetries of the rate ^4cp an d the angular coefficients [7] 

A < s WTf§? for ' = 3 ' 6 ' 9 ' A - s WTwB for ! = 4 ' 5 ' 7 ' 8 (4) 

based on naive factorisation was extended using QCDF to the low-g 2 region by taking into 
account the NLO QCD corrections for the strong phase 5 S [10]. Both A^g vanish in the 
SM at leading order in QCDF since they are proportional to the interference of SM and 
chirality flipped operators, making them ideal probes of the latter. Furthermore, A^g are 
subject to cancellations at leading order in QCD, such that the NLO QCD corrections give 
large corrections. The (^-integrated SM predictions of the CP asymmetries are summarised in 
table 1. They are, as expected, tiny and below a percent. A model-independent analysis in the 
presence of the BSM contributions to the SM and the chirality flipped operators 0' 7 9 10 shows 
substantial room for enhancement. The T-odd asymmetries A^g 9 can be of order one. However, 
the experimental uncertainties at LHCb are large for CP-asymmetries, even when considering 
improved normalisations as discussed for Aq 2s and A% [27]. 

The rate-normalised CP-averaged quantities Sf ~ (J/ + J[) constitute another important 
set of observables, of which the CP-odd ones (i = 5, 6, 8, 9) require i?-tagged samples. The 
Si have been analysed in the SM, model-independently and model-dependently in great detail 
[12, 13]. Especially, the lepton forward-backward asymmetry ^4fb (related to Sg) and S§ could 
be measured rather precisely with early LHCb data of 2 fb _1 
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Figure 2. Comparison of BaBar [yellow], 
Belle [red] and CDF [black] data points with 
SM predictions (blue) for B% K*° + It 
distributions [4]: dB/dq 2 , Afb and Fl. The 
vertical grey bands are vetoed in experiment. 



It should be noted that [26, 27] estimate uncertainties due to unknown subleading 
0(Aqcd/Ek*) contributions to both form factors and B — > K* + ££ amplitudes by power 
counting. Contrary, [12] does not make use of the large recoil symmetry relations of the form 
factors in the LO QCD contribution of the said amplitudes. Instead, all seven QCD form factors 
determined from light cone sum rules (LCSR) [29] are used, thus accounting for the associated 
unknown subleading 0(Aqqd/Ek*) contributions. The direct application of the LCSR results 
reduces the form factor uncertainties [12] due to more efficient cancellations. 

Recently, soft-gluon emission effects due to (cc)-resonances were found to extend down to the 
low-q 2 region [30] affecting the ^-differential rate at most about +15 % at q 2 = 6 GeV 2 . Due 
to this relatively large shift such effects should be investigated further and should be considered 
in future studies of the exclusive rare decays. 



2.2. High-q 2 / Low recoil region 

Assuming quark-hadron duality, the low recoil region can be studied using an OPE in 
combination with HQET to treat the long-distance contributions from current-current and QCD- 
penguin operators [19, 20]. The authors of [20] perform the OPE without subsequent matching to 
HQET contrary to [19], among further technical details. Duality violating corrections to the OPE 
are estimated in [20] and found to be small. In [19] it is shown, that at the leading order in the 
Aqcd/Q, Q = {Vq 2 , mb} expansion 8 all hadronic matrix elements can be expressed through the 
3 QCD form factors V,A\ t 2- As a consequence, a universal short-distance dependence emerges 
for all transversity amplitudes (a = 0, _L, ||) [4] 

A L a ' R oc c L , R x f a , c L , R = (Cf t Cio) + it^Cf, (5) 

with the (effective) short-distance couplings Cf*| 10 and the form factor terms f a (a = 0, _L, ||) 
which are linear in V,A±^. The factor k ~ 1 is known including NLO in QCD from matching 

8 Note, that the subleading order is additionally suppressed by a s . 



QCD onto HQET [19]. The term oc Cj S is subject to uncertainties of order AqcdA^6 due to the 
HQET form factor relations, which are additionally suppressed in the SM by \C-jjCg\ ~ 0.1. As 
a consequence of (5) J7 = Js = J9 = and the remaining angular observables 

Jf ~ U x = 2 Pl fl (2Ji + J 3 )~U 2 = 2 Pl fl, (2JI-J 3 )~U 3 = 2 Pl fl (6) 
Ji ~ U4 = 2pi/ /||, J5 ~ U 5 = 4/? 2 /o/_l, J| ~ ^6 = 4p 2 /||/±, 

depend only on the two short-distance (dominated) combinations p\ = (|cr| 2 + |cl| 2 )/2 and 
P2 = (|cr| 2 — |cl| 2 )/4 [4]. From (6) one obtains the three long-distance free ratios 



r (i) = Ui = 1 H (2) = U5 =2 P2 R (3) = U% =2 P2 

T VU^Ui ' T VU^Th Pi T VU^Th Pi 



After integration over q 2 G [14.0, 19.2] GeV 2 one obtains in the SM 

(fl-W) = +0.997 ± 0.003, (H^) = -0.972 ± 0.010, {H^} = -0.958 ± 0.010. (8) 

Details concerning the definition of (. . .) and the theoretical uncertainties are given in [4] as well 
as the results for B,A F b,Fl and A^' 3 '^. We stress that (H^) rj 1 holds model-independently. 
Hence, deviations from this prediction test the validity of the OPE framework. Furthermore, 
the following short-distance free ratios 



fo fTh U\ U 4 U 5 f fTh fjL \U~2 ^ 



/11 \U 3 U A U 3 U 6 ' /j_ Vt/ 2 ' f\\ \u 3 U A 



(9) 



allow to test the B — > V form factors at high-q 2 , such as obtained from lattice QCD, using 
experimental data [4]. 

A comparison of the SM predictions and the available data from BaBar, Belle and CDF for 
low- and high-g 2 regions is shown in figure 2 where the main uncertainties (shaded bands) arise 
from the form factors and the missing subleading corrections of order Aqcd / mb as estimated 
by power counting [4]. The Belle and CDF measurements of the three low- and high-(7 2 bins 
q 2 £ [1,6], [14.18, 16.0], [16.0, 19.21] GeV 2 of B and A FB (and q 2 G [1,6] GeV 2 of F L ) yield the 
first constraints on the short-distance couplings Cg t w of the SM operators Og,io beyond naive 
factorisation. The constraints for real-valued Cg,io are shown in figure 3 fixing C7 = +Cf M 
(results for C-j = — Cf M are available in [4]) when employing low-g 2 data only and when adding 
high-g 2 data 9 , respectively. Figure 4 shows |Cio| vs. \Cg\ for complex-valued Cg : \o and results of 
(^9 vs. |Cg|), (4>io vs. |Cio|) and (</>iq vs. 0g) can be found in [5]. Both figures demonstrate the 
strong constraining power of the high-g 2 region on Cg 5 io, and consistency with low-g 2 data. 

CP asymmetries with improved short-distance sensitivity can be designed in analogy to 
Ht' 2 ' 3 \ since j{ ~ p a with p a = p a [Sw — > —8w]- It is convenient to consider the three 
CP-asymmetries 

P2 _ fh 

(1) _ Pi ~ Pi (2) _ pi ft (3) _ „ P2 ~ P2 nm 

a cp- pi+ ^ a cp-p^ + |> a cp- z pi + ^ l iU J 



Using an extrapolation of the B — > V form factors from low-g 2 QCD LCSR [29] . 
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Figure 3. Constraints on 
real- valued £9,10 at 68%, 
95% and 99% CL from low- 
q 2 bins of B° -> K*°££ and 
B -> X S U [left] and when 
adding high-g 2 bins [right] 
[4, 5]. Here C 7 = +Cf M is 
fixed. The square marks the 
SM value of (C 9 ,C W ). 



Figure 4. Constraints 
on |C 9i io| at 68% and 95% 
CL for complex-valued £9,10 
from low-g 2 bins of B^ — > 
K*°l£ and B -> X S U 
[left] and when adding high- 
ly 2 bins [right] (see [5] for 
details). The square marks 
the SM value of (C 9 ,Ci ). 



with various possibilities to measure them [5] . Moreover, a C p can be extracted from an untagged 
i?-meson sample. There are no T-odd CP-asymmetries at low- recoil due to (6): J7,s,9 = 0. In 
the SM |aQp 2,3 ^|sM ^ 10~ 4 whereas a model-independent analysis for complex-valued Cg^o allows 

(1 3) (2) 

values up to \(a C p )| < 0.2 and [(a^p)! ^ s presently unconstrained (see [5] for details). 
3. CP asymmetries without tagging 

The decays B S ,B S — > <f>{—> K~ K + ) + U require to account for the mixing of the initial B®- 
and i?^-mesons before decaying into the common final state. Consequently, their observables 
depend additionally on the B s — B s mixing phase <3? s and the lifetime difference Ar s , which are 
currently being analysed at the Tevatron, for instance in B s , B s — > </>(— )■ K + K~) + J > U). 
In this respect, the CP-odd property of the observables J\ with i = 5, 6, 8, 9 allows to measure 
the associated CP-asymmetries from an untagged and time-integrated data set. 

The CP-asymmetries A™ lx for i = 6, 9 and Af mix for i = 5,8 have been worked out at 
low-g 2 [10]. They match the corresponding Ai and Af of the flavour-specific unmixed decays 
for Ar s — > 0. There the CP-averaged rate has been used as normalisation. An improved 
cancellation of the hadronic uncertainties might be achieved along the lines of the unmixed case 
by choosing the normalisations proposed for Aq 2s and A% [27]. 

At high-q 2 the CP-asymmetries associated with ■/•?, i = 8, 9 vanish and those with i = 5, 6 



can be related to the same short-distance dominated a™™ using the normalisations as indicated 
in reference [5]. Notably, a™™ is independent of the sign of Ar s and the dependence on & s is 
rather weak for small AT S /F S . Experimentally, Ar s /r s ~ O(0. 1). 

4. Conclusion 

The angular analysis of the exclusive B° d ->• K~n + ) +U and B S ,B S ^ <£(->• K-R + ) + U 

decays with 4-body final states offers a rich phenomenology to test the short-distance flavour 
couplings in a large number of observables. The LHCb experiment is expected to improve the 
experimental situation tremendously in the near future, but also the upcoming final analysis of 
the complete BaBar, Belle and CDF data sets will contribute to the exploration of the borders 
of the SM. 

The low- and high-g 2 regions are accessible to power expansions which reveal symmetries of 
QCD and indicate suitable combinations of observables to reduce the hadronic uncertainties. 
This allows for precise theory predictions for the exclusive decays. 

The low-g 2 region is theoretically (QCDF, SCET) well understood and effects of (cc)- 
resonances can be estimated. It offers many interesting observables which wait for confrontation 
with data. 

The high-g 2 region is based on an OPE relying on quark-hadron duality. Violations of the 
latter could be signaled experimentally by deviations from = 1. Currently, the B — > V 
form factors used at high-g 2 are extrapolated from the low-q 2 region using a dipole formula, but 
lattice calculations for the high-g 2 regions are in progress [31]. There are two "long-distance 

(2 3) 

free" ratios Hjp ' which provide tests of the SM and, moreover, several "short-distance free" 
ratios which offer direct comparison of the lattice results for ratios of form factors with data. 

A new flavour tool "EOS" aimed at the evaluation of the observables covered in this article 
is developed at TU Dortmund [32] whose first stable release is expected during 2011. 
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